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Abstract. Origin and properties of non-Lorentzian spectral lines in linear chains 
of nanospheres are discussed. The lines are shown to be super-exponentially narrow 
with the characteristic width oc exp[— C(h/a) 3 } where C is a numerical constant, h 
the spacing between the nanospheres in the chain and a the sphere radius. The fine 
structure of these spectral lines is also investigated. 
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One-dimensional periodic chains (ODPC) of metallic nanospheres have attracted 
significant recent attention due to their unusual optical properties. Although the general 
theoretical framework for analyzing electromagnetic interactions in ODPC has been built 
a decade ago pQ, the recent dramatic advances in nanofabrication have reinvigorated the 
interest in ODPC, which, in turn, has led to several new results of high experimental 
relevancy. In particular, radiatively non-decaying surface plasmons (SPs) in ODPC with 
possible applications to building novel lasers were discussed in Ref. [2|; unusual shifts of 
plasmon resonance frequencies were found in Ref. [Sj and a dramatic narrowing of SP 
spectral lines was found in Ref. jH |3] in finite chains of moderate length. In this letter I 
show that two of these phenomena (unusual shifts and narrowing of SP spectral lines) are 
directly related to a logarithmic divergence of dipole sums (electromagnetic eigenvalues) 
- a theoretical interpretation that has not been given so far. The SPs that can be 
excited as a result of this divergence posses highly unusual properties. In particular, the 
resonance line-shapes are essentially non-Lorentzian and are characterized by a vanishing 
integral weight. This is in a sharp contrast to spectral line broadening or narrowing due 
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Figure 1. Sketch of the ODPC excitation by a linearly polarized plane wave. 



to change in the decay rate. Another consequence which has not been previously noticed 
is that each narrow resonance is paired with an even more narrow spectral hole. An 
interesting new property discussed below is that the narrow collective SP resonances 
can be excited in ODPCs even when the distance between neighboring spheres is much 
larger than the sphere diameters. However, the resonances become so narrow in this 
case that, unless one takes special care, it is extremely unlikely to notice them in any 
numerical or experimental investigation. 

We start with a brief summary of the underlying physics. The geometry of a 
ODPC excited by a plane wave of the form E exp(zk • r) is illustrated in Fig. ^ 
Here the polarization vector Eo is perpendicular to the chain, the incidence angle is 
denoted by 9, lattice spacing by h and the sphere radius by a. Each sphere is assumed 
to be polarizable and characterized by the dipole polarizability a. We work in the 
approximation introduced by Doyle jH] in which each sphere is treated as an elementary 
dipole located at its center but is characterized by non-quasistatic polarizability a which 
is calculated from the coefficient a\ of the Mie theory [2j : 



where ipi(x) and are the Riccati-Bessel functions, m = y/e is the complex refractive 
index of the spheres and k = uj/c = 2n / A is the wave number of the incident wave. The 
above approximation allows one to include spheres which are not small compared to 
the incident wavelength A while staying within the purely dipole theory. The higher 
multipole interactions of the spheres, as well as the input of higher multipoles to the 
optical cross sections, are ignored in this approximation. Note that the polarizability 
a defined by (JTJ) is the exact dipole polarizability with respect to excitation by a plane 
wave, but not by secondary waves scattered by spheres in the chain. However, the 
dipole approximation was shown to be very accurate when h is the order of or larger 
than 2a (which is the case discussed below) by direct comparison with a converged 
T- matrix solution 4J. In general, it is known that short-range multipole interactions 
of orders higher than the first (dipole) do not play a significant role for transverse 
electromagnetic excitations of finite or infinite linear arrays of interacting spheres even 



a = 



3i m^iimka^'^ka) — tpi(ka)ip[(mka) 



(1) 



2k 3 mi/ji(mka)C,[(ka) — ^i(ka)ip' 1 (mka) 
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Figure 2. (a) Function ReS'(x) for 6 — 0. (b) Re(a 3 /a) calculated for silver for 
different values of a; the curves a = 5nm and a = lOnm correspond to the quasi-static 
limit and are indistinguishable. 



when the spheres are in close proximity [HI Ej . Physically, one can argue that the short- 
range interaction is not important for transverse excitations because it does not lead to 
an electric current along the chain (in a sharp contrast to the longitudinal excitations 
which are not discussed in the letter). 

In the approximation formulated above, each sphere is characterized by a dipole 
moment which, in the case of geometry shown in Fig. Q is collinear with the polarization 
vector Eo and has the amplitude d n , n being the index which labels spheres in the chain. 
The amplitudes d n are coupled to the external field and to each other by the coupled- 
dipole equation [l] 



d n = a 



E exp (ikn sin 9) + 2J W n - n '(kh)d n i 



Here W n _ n / is the dipole interaction term given by 



W n (x) = k s 



' —1 • —2 — 3 " 

X IX X 



+ 



exp(2x|n|) . 



\ \ n \ \ n \ \ n \ 
The coupled-dipole equation (J2J is easily solved to yield 
a 3 Eo exp(ikhn sin 9) 



d 



a 3 1 a — (ka) 3 S(kh) 



(2) 



(3) 



(4) 
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where the dimensionless quantity (ka) 3 S(kh) is the dipole sum that characterizes 
excitation of the SP with the wave number q = ksin8. The function S(x) is given 
by 

S(x) = 2 f h - — — - — - — rr | exp(mx) cos(nx sin 9) . (5) 

~^ \xn (xn) 2 (xny ) 

It can be seen that the above series diverges logarithmically when (1 ± sm9)kh = 2nl, 
I being an integer. It is convenient to separate the sum into two parts: S(x) = 
Si(x) + S 2 (x) where Si is given by 

n i \ o cos(ra) cosfnxsin^) 1 r . 2 -i 

Si(z) = 2V — - = In 4(cosx - cos(xsm6>)) (6) 

and diverges when cosx = cos(xsin#) while S^x) is the reminder of series (0) and 
converges for all values of parameters. For simplicity, we will assume everywhere below 
normal incidence (6 = 0) which was also the case considered in Ref. jlj. Then the 
divergence takes place when h = XI. 

The specific extinction cross section per one sphere is given by pQ 

Anka 3 

° e m a 3 /a-{ka) s S{kh) ' 1 ' 

Optical resonances occur when the real part of the denominator in the above expression 
vanishes. In Fig.[^a) we plot ReS'(a;) for 6 = 0. The sharp peaks in the plot correspond 
to the points where ReS" diverges. In Fig. |2Jb) we also plot Re(a 3 /a) calculated 
according to for silver. Interpolated experimental dielectric function from was 
used in calculations. It can be seen that in an isolated sphere, the SP (Frohlich) 
resonance takes place in the interval 350nm < A < 380nm, depending on the value 
of a. The resonant wavelength is obtained from Re(a 3 /a) = 0. Above the Frohlich 
resonance (at smaller wavelengths), the spectral variable Re(a 3 /a) becomes negative. 
(Here we ignore the region A < 320nm where no resonance excitation can take place 
due to the strong interband absorption.) Therefore, in order to excite a SP in an 
interacting ODPC in this spectral region, the variable ReS* must be also negative. As 
can be seen from Fig. Efa), this happens for sufficiently small values of the parameter 
x = kh and corresponds to the conventional blue shift of the transverse electromagnetic 
oscillations, which is well known and can be described by the quasistatic interaction 
However, below the Frohlich resonance (at larger wavelength), the spectral parameter 
Re(a 3 /a) is positive. Therefore, in order to excite an SP in this spectral region, the 
variable ReS* must be also positive. Obviously, this requirement is fulfilled near the 
points of divergence of ReS*. Thus, if h/X is close to an integer, the transverse collective 
oscillations of the chain are shifted to the red from the Frohlich wavelength, contrary 
to the usual case. Quite remarkably, the collective resonance can take place even if 
a << h and (ka) 3 << 1. Indeed, no matter how small ka is, the resonant condition can 
always be satisfied sufficiently close to the point A = h/l. Below, we focus on the first 
of these resonances, which corresponds to A ~ h and can be experimentally observed in 
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Figure 3. (a) Extinction efiiciency Q e for ft, = 500nm and a — 50nm. (b) Illustration 
of cancellation of the real part of the denominator in (JJJ . (c) Spectral dependence of 
the imaginary part of the denominator, (d) The narrow resonance shown in the panel 
(a) fully resolved. 



metal ODPCs in the visible and IR spectral regions below the Frohlich frequency of an 
isolated sphere. 

In Fig. Of a) we plot the dimensionless extinction efficiency Q e = a e /A7rka 3 for 
the following values of parameters: a = 50nm and h = 500nm. The sharp resonance 
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corresponding to A ~ h is clearly visible. The curve is very close to the one shown 
in Fig. 1(a) of Ref. [3]. In Fig. 0(b) the origin of the sharp resonance is illustrated. 
Namely, we see that in a very narrow spectral interval near A = 500nm real parts 
of the two terms in the denominator of (JJJ) cancel each other. From the analysis 
of Fig. 0(b), it is obvious that there should be, in fact, two closely spaced narrow 
resonances separated by a spectral hole. However, the first of the two resonances is 
suppressed due to radiative losses. The imaginary part of the denominator of (JJJ) 
is plotted in Fig. 0(c)- As was pointed out in Ref. [T|, the imaginary part of S(x) 
experiences a jump at x = 2tt. However, no exact cancellation of the imaginary part 
of the denominator can take place. For the geometry considered here, it can be shown 
that lm[(ka) 3 S (kh)] > — 5(ka) 3 /12 while Im(a 3 /a) < — 2(fca) 3 /3, so that the imaginary 
part of the denominator is < — (ka) 3 /4, with the equality taking place for nonabsorbing 
materials with Ime = 0. Finally, In Fig. 0(d) the sharp resonance seen in Fig. 0(a) is 
completely resolved. The narrow spectral hole located exactly at A = h can be also seen 
in this figure. 

Let us estimate the width and amplitude of the narrow resonances occurring due 
to the divergence of S(x). We define the width of the resonance as the distance 
from the resonance wavelength A r , determined from the condition Re[a 3 /o;(A r ) — 
(2iia/ \ r ) 3 S (27th/ \ r )} = to the center of the spectral hole at A = h] thus, AA = |A r — h\. 
We estimate A r assuming that the dominant contribution to KeS comes from the 
logarithmically diverging term (JHJ). We also assume that A A << 2n and expand the 
argument of the cosine in (JHJ) near the point x = kh = 2tt, which leads to the following 
estimate: 



where C = Re[a 3 /a(A = h)} is a constant of the order of unity. For example, for 
h = 500nm, C ~ 0.5. Using a = 50nm, we obtain from (|SJ) AA ~ 0.14nm in agreement 
with Fig. 0(d). Thus, the width of the resonance is completely determined by the 
geometrical factors (the ratio h/a) and is not in any way controlled by relaxation. The 
latter, however, influences amplitude of the resonance. Indeed, the maximum value of 
Q e in the peak is given by l/Im[a 3 /a; — (ka) 3 S(kh)). For the geometry considered here, 
it can be verified that this value can not be greater than (/i/Va) 2 , which is the limit 
for nonabsorbing material. However, for strongly absorbing materials amplitude of the 
resonance can become negligibly small. 

Since the amplitude of the narrow resonances does not increase when the width 
decreases super-exponentially, it is impossible to effectively excite these resonances by 
a near-field probe. For example, consider the case when a single sphere (say, n = 0) 
is excited by a near-field microscope tip of small aperture. Then the coupled-dipole 
equation for the amplitudes d n can be solved by Fourier transformation: 




(8) 




(9) 
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where S(kh,qh) is given by JHJ) in which xsinO in the argument of cosine must 
be formally substituted by qh and the variable x in the reminder of the formula 
substituted by kh. Function S(kh, qh) diverges logarithmically when cos(kh) = cos(qh); 
in particular, if kh = 2tt as in the examples considered above, the only point of 
divergence within the integration interval is q = 0. Similarly to resonances in extinction 
spectra, this resonance is super-exponentially narrow in the SP wavenumber q, and its 
input into the above integral is negligible. We emphasize again that the resonances 
discussed here are essentially non-Lorentzian, and the conditions for applicability of the 
quasiparticle pole approximation, which under normal circumstances would properly 
describe coupling of the near-field probe to SPs, are severely violated. 

The narrow resonance in Fig. Efd) was obtained for h = 500nm and a = 50nm. We 
now show that the narrow resonances also exist for smaller values of a and larger ratios 
h/a. To this end, we plot the extinction efficiency Q e for h = 500nm and a = 45nm 
[Fig Ufa)] and a = 40nm [Fig 13(b)] . The dielectric function does not vary noticeably 
over the narrow spectral range shown in Fig. 0] and was therefore taken to be constant, 
e = —8.5 + 0.762, which corresponds to the experimental value at A = h — 500nm given 
in Ref. [10J. The narrow non-Lorentzian resonances are well manifested in Fig. |3J The 
central spectral hole is resolved to some degree in Fig. 0(a) but is shown only as a vertical 
line in Fig. Efb). Obviously, it is impossible to resolve the spectral holes completely, 
since they do not have the Lorentzian structure and are non-differentiable at the point 
A = h. Note that the resonances shown in this figure are much more symmetrical with 
respect to the point A = 500nm than the one shown in Fig. 13(d). This is due to the 
fact that with increasing the ratio h/a, the influence of radiative losses on the shape of 
resonance lines decreases. For even larger values of h/a, the resonances quickly become 
extremely narrow but do not disappear completely, at least in chains of sufficient length. 

It is interesting to consider the possibility of narrow resonances in situations when 
the dipole approximation is not applicable, i.e., for spheres in close proximity. It can 
be shown that the resonances discussed in this letter do not disappear or get broadened 
when the full multipole interaction is taken into account. Furthermore, the spherical 
shape of the particles is also not fundamental because the phenomenon discussed here 
originates due to long-range interaction in ODPC while the higher multipole interaction 
is short-range. This conclusion is in agreement with the numerical study of extinction 
spectra of periodic chains of cylindrical disks 5J which were shown to have sharp 
resonances similar to those found in chains of spheres. Generalization to two-dimensional 
arrays of particles is also possible. Two-dimensional periodically modulated structures 
have also attracted significant recent attention, with the possible application including 
random lasers [2|, development of novel chemical and biological sensors [3] and the study 
of anomalous optical transmission through metal films [T2] . 

To conclude this letter, we discuss several factors that contribute to broadening 
of the spectral lines discussed above. The most important factor is the finite length 
of a chain, since the divergence of the dipole sums is logarithmic. As was mentioned 
above, narrow resonances very close to those in infinite chains were obtained in Ref. [3] 
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Figure 4. (a) Narrow resonances for h = 500nm and a — 45nm (a) and a = 40nm 
(b). 



for h/a = 10 and only 50 spheres; however, observing more narrow resonances with 
h/a > 10 will require a substantially larger number of spheres. One possible solution 
to this problem is to place the ODPC into a circular optical fiber. Disorder is another 
important factor. Numerical simulations in finite chains (400 particles) [Sj revealed that 
random uncorrelated displacements of particles with the amplitude of ~ O.lh do not 
noticeably change the resonance lineshape. This is an expected result for short-range 
disorder, i.e., the disorder with the correlation length of one or few lattice spacings. 
However, disorder with long-range correlations can result is much stronger changes in 
the resonance lineshapes. Further, the account of nonlocality of the dielectric response 
will not alter the nature of positive interference (synchronism) which results in the 
logarithmic divergences, and is not expected to broaden the narrow spectral lines. 
The two physical phenomena whose effects on the spectral lines discussed here are 
difficult to predict are the nonlinearity of the optical response (e.g., Kerr-type third- 
order nonlinearity) and quantum effects. These effects will be the subject of future 
work. 
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